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Abstract.  We  report  on  ~  120  days  of  observations  of 
Cygnus  X-1  with  OSSE  onboard  the  Compton  Observa¬ 
tory.  Emission  is  detected  in  the  range  50  keV  to  1  MeV, 
and  we  find  evidence  for  a  continuum  of  hard  X-ray  spec¬ 
tra  over  a  broad  range  of  fluxes  rather  than  the  existence 
of  distinct  flux  states.  Comparisons  of  the  source  spec¬ 
tra  with  various  theoretical  models  show  that  an  expo¬ 
nentially  truncated  power  law  best  describes  the  average 
spectrum  in  the  OSSE  energy  band.  Although  we  have 
measured  a  new  minimum  in  the  hard  X-ray  flux  from  the 
source,  no  evidence  was  found  for  either  a  broad  1  MeV 
feature  or  a  narrow  511  keV  line  previously  reported  in 
association  with  a  low  flux  state.  Upper  limits  on  such 
emission  features  are  an  order  of  magnitude  lower  than 
earlier  reported  detections.  The  5.6-day  periodicity  of  the 
source  measured  at  optical  wavelengths  was  not  detected 
with  a  sensitivity  to  the  rms  modulation  fraction  of  5%  in 
the  60-140  keV  energy  band. 

Key  words:  Gamma  Rays:Observations-  Stars:  Individ¬ 
ual:  Cyg  X-1  -  Black  hole  physics 


1.  Introduction 

Cyg  X-1,  one  of  the  brightest  and  most  extensively  studied 
gamma  ray  sources  in  the  sky,  is  a  5.6-day  binary  system 
consisting  of  a  blue  supergiant  (HDE  226868)  and  a  com¬ 
pact  companion  of  mass  ~6.3  Mq  (Gies  et  al.  1986,  Dolan 
et  al.  1989).  Since  the  compact  object  has  a  mass  greater 
than  3  Mq,  Gyg  X-1  is  generally  considered  to  be  a  prime 
black  hole  candidate,  and  it  has  become  the  prototype 
for  back  hole  characterization  and  modeling.  The  possible 
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presence  of  features  in  the  spectrum  of  Gyg  X-1  has  been 
suggested  as  a  useful  criterion  for  distinguishing  black  hole 
systems  from  those  containing  neutron  stars  (Liang  et  al. 
1992).  In  particular,  features  such  as  soft  X-ray  excesses  or 
gamma-ray  bumps  in  the  hard  X-ray /gamma-ray  energy 
band  have  been  proposed  as  signatures  of  black  holes,  giv¬ 
ing  special  importance  to  an  improved  knowledge  of  the 
detailed  shape  of  the  Gyg  X-1  spectrum  at  high  energies. 

The  high-energy  radiation  from  Gygnus  X-1  is  thought 
to  be  emission  from  an  accretion  disk  surrounding  a  black 
hole  (Shapiro,  Lightman  &:  Eardly  1976).  The  soft  1-10 
keV  X-ray  spectrum  has  traditionally  been  described  by  a 
dominant  low  state  with  occasional  periods  of  high  state 
emission  (e.g.,  Priedhorsky  et  al.  1983).  The  hard  X-ray 
flux  is  well  represented  by  a  thermal  Gomptonization  spec¬ 
trum  (Sunyaev  and  Titarchuk  1980).  Using  HEAO-3  data. 
Ling  et  al.  (1987)  identified  three  distinct  levels  for  the  45- 
140  keV  flux  occuring  during  the  low  X-ray  state,  which 
they  labeled  7I,  72  and  78,  with  7I  being  the  state  with 
the  lowest  45-140  keV  flux. A  gamma-ray  excess,  charac¬ 
terized  by  a  1.2  MeV  wide  gaussian  at  1  MeV  with  a  line 
flux  of  1.6  X  10“^  ph  cm“^  s“^,  was  detected  in  fall  1979 
when  Gyg  X-1  was  in  the  7I  state  (Ling  et  al.  1987).  A 
weak  2(7  detection  of  a  narrow  511  peak  was  later  reported 
from  analysis  of  the  HEAO  data  during  the  same  episode 
(Ling  et  al.  1989).  Other  groups  have  also  reported  mea¬ 
surements  of  MeV  emission  from  Gyg  X-1  (see  Owens  and 
McGonnell  1992  for  a  review  of  the  MeV  observations). 
Recently,  Harris  et  al.  (1993)  derived  upper  limits  for  flar¬ 
ing  episodes  of  12  or  more  days  from  the  long-term  data 
base  of  the  gamma  ray  spectrometer  on  SMM,  and  Mc¬ 
Gonnell  et  al.  (1994)  reported  no  evidence  for  MeV  fea¬ 
tures  from  the  Gyg  X-1  in  data  from  GOMPTEL. 

In  this  paper,  we  present  results  of  measurements  of 
Gyg  X-1  using  OSSE.  Given  the  importance  of  testing 
the  models  of  the  source  and  resolving  the  conflicting  na- 
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ture  of  some  of  the  previous  measurements,  OSSE  devoted 
122  days  to  observing  the  source  at  different  intensity  lev¬ 
els  during  the  first  four  years  of  the  CGRO  mission.  The 
shape  of  the  spectrum,  the  search  for  narrow  and  broad 
spectral  features,  the  light  curve  of  daily  flux  averages, 
the  flux  levels  of  the  source,  and  correlations  of  hardness 
ratios  are  addressed. 


2.  Observations 

The  OSSE  experiment  observed  Cyg  X-1  during  17  inter¬ 
vals,  each  spanning  one  or  more  days  in  the  first  four  years 
of  the  CGRO  mission,  for  a  total  of  ~2.4  x  10®  live  seconds 
of  background-subtracted  spectroscopy  data.  The  observa¬ 
tions  were  during  the  following  CGRO  viewing  periods:  2, 
7,  15,  203,  209,  212,  221,  223,  303,  318,  328,  331,  331.5, 
332,  333,  419.5,  420.  These  data  are  not  a  random  sample 
of  the  flux  levels  of  the  source  because  four  of  the  obser¬ 
vations,  viewing  periods  212,  221,  223  and  318,  were  in  re¬ 
sponse  to  a  low  hard  X-ray  state  detected  by  the  Burst  and 
Transient  Source  Experiment  (BATSE).  We  performed  a 
standard  spectroscopic  analysis  of  the  data  (Johnson  et 
al.  1993)  for  each  detector  and  each  day  separately,  and 
the  spectra  were  added  at  the  end  of  the  analysis.  The 
timing  analysis  of  these  data  at  mHz  to  kHz  frequencies  is 
not  considered  here  and  will  be  presented  elsewhere  (see 
Bridgman  et  al.  1994  for  an  early  analysis). 
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Fig.  1.  Spectrum  of  the  sum  of  all  Cyg  X-1  data  with  the 
best  fit  exponentially  trnncated  power  law  (sohd)  and  thermal 
Comptonization  model  (dotted). 


3.  Temporal  Analysis 

The  integrated  flux  in  the  45-140  keV  band  was  monitored 
for  each  day  of  data  and  a  continuum  of  intensities  rather 
than  obvious  discrete  states  was  detected.  This  becomes 


clearer  when  the  daily  flux  is  observed  over  a  long  pe¬ 
riod  of  time,  as  is  done  by  BATSE  (Paciesas  et  al.  1994). 
Compared  with  the  historic  7I,  72  and  73  levels,  the  ob¬ 
servations  of  February  1994  (CGRO  viewing  period  318) 
revealed  a  flux  intensity  of  ~  1/5  the  previously  reported 
lowest  flux  in  this  band.  The  source  has  therefore  a  greater 
range  of  intensities  than  previously  measured.  Note,  how¬ 
ever,  that  no  OSSE  measurements  have  been  performed 
at  intensities  comparable  to  the  78  level. 

The  light  curve  shown  can  be  searched  for  periodici¬ 
ties,  in  particular  for  the  5.6-day  orbital  periodicity.  The 
data  set  is  sparse,  with  only  ~100  days  of  data  out  of 
a  1024-day  interval,  with  the  contiguous  daily  data  sets 
varying  from  two  days  to  two  weeks  in  length.  We  calcu¬ 
lated  the  Fast  Fourier  Transform  for  the  first  1024  days 
of  the  mission  after  we  subtracted  a  fifth  order  polyno¬ 
mial  in  time  to  detrend  the  data. We  find  no  evidence  for 
periodicity  in  this  sparse  data  set,  other  than  a  peak  at 
very  low  frequency  due  to  the  data  sampling.  The  upper 
limit  to  the  amplitude  of  a  5.6-day  periodicity  is  0.005 
7cm“^  s“^  in  the  60-140  keV  band,  corresponding  to  an 
rms  modulation  fraction  of  <  5%. 


4.  Spectral  Analysis 

The  spectrum  for  the  sum  of  all  data  through  May  1995 
is  shown  in  Figure  1.  The  spectrum  is  smooth,  with¬ 
out  breaks  or  line  features.  It  extends  from  60  keV, 
the  lower  threshold  used  in  this  spectral  analysis,  to  ~1 
MeV,  where  the  flux  level  falls  below  the  OSSE  sen¬ 
sitivity.  The  high  statistical  precision  of  the  Cyg  X-1 
spectrum  permits  differentiation  among  theoretical  mod¬ 
els  to  much  higher  energies  than  previously  possible.  We 
use  the  four  different  functional  forms  to  fit  the  data: 
a  two-parameter  thermal  bremsstrahlung  model  without 
Gaunt  factor(flux  at  100  keV,  =0.476  7cm“^s“^MeV“^, 
kT=108  keV),  the  three-parameter  thermal  Comptoniza¬ 
tion  spectrum  of  Sunyaev  and  Titarchuk  (1980)(flux  at 
100  keV=  0.465  7cm“^s“^MeV“^,  kT=65  keV,  r=1.94),  a 
three-parameter  exponentially  truncated  power  law  (flux 
at  100  keV=  0.470  7cm“^s“^MeV“^,  Ecutoff  =  ^^^  keV, 
photon  index=1.39),  and  we  also  include  the  possibility 
that  some  fraction  of  the  emission  described  by  an  ex- 
pentially  truncated  power  law  is  incident  on  a  cold  opti¬ 
cally  thick  scattering  medium  (flux  at  100  keV=  0.473  ph 
cm“^s“^MeV“^,  Ecutoff  =  ^^^  keV,  photon  index=0.96, 
reflection=1.5).  The  effects  of  Compton  reflection  are  con¬ 
tained  in  a  code  provided  by  Zdziarski  (1995)  to  perform 
the  analysis.  The  model  contains  a  reflection  parameter 
that  is  proportional  to  the  solid  angle  covered  by  the  re¬ 
flection  medium.  All  fits  were  over  the  energy  range  from 
60  keV  to  4  MeV. 

The  model  with  an  exponentially  truncated  power  law 
gives  the  best  fit  to  the  data.  The  thermal  bremsstrahlung 
and  the  thermal  Comptonization  models  both  diverge 
from  the  measured  spectrum  at  higher  energies  where  they 
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underestimate  the  flux  from  the  source.  The  exponentially 
truncated  power  law  with  reflection  produces  a  slightly 
worse  fit  to  the  data,  and  the  fit  improves  as  the  amount 
of  reflection  decreases.  The  difference  in  the  photon  spec¬ 
tra  between  the  model  with  and  without  reflection  oc¬ 
curs  mostly  at  the  lower  end  of  the  OSSE  energy  range, 
where  the  systematic  errors  in  the  instrument  efficiency 
and  calibration  are  largest.  The  confidence  in  rejecting 
models  based  solely  on  the  lowest  energy  bins  is  there¬ 
fore  smaller  than  the  confidence  in  rejecting  models  which 
diverge  at  higher  energies.  The  exponentially-truncated 
power  law  is  also  the  best  model  when  each  viewing  period 
is  summed  and  fitted  separately.  The  optically  thin  ther- 
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Fig.  2.  Temperature  versus  45-140  keV  integrated  flux,  for 
each  day,  obtained  using  an  optically  thin  bremsstrahlung 
model. 

mal  bremsstrahlung  model  still  provides  a  good  fit  to  the 
daily  spectra.  This  simple  two-parameter  model  is  used  to 
study  the  temperature-intensity  dependence  of  the  source. 
Figure  2  shows  that  the  temperature  of  the  source  and  the 
intensity  vary  within  a  limited  range,  except  for  a  few  ex¬ 
cursions  to  the  low  temperature-  low  amplitude  regime. 
During  the  very  low  flux  observation  of  February  1994, 
the  temperature  of  the  source  was  ~  70  keV  compared 
with  110  keV  for  the  other  observations.  Note  that  the 
scatter  in  the  points  on  this  figure  is  much  larger  than  the 
statistical  or  systematic  errors.  Alternatively  fitting  the 
daily  data  with  the  thermal  Comptonization  model  de¬ 
rived  by  Sunyaev  and  Titarchuck  (1980),  we  find  a  similar 
scatter  in  the  optical  depth  parameter,  with  the  tempera¬ 
ture  fairly  constant  with  a  mean  of  65  keV  and  a  standard 
deviation  of  less  than  5  keV. 

In  Figure  3,  we  display  spectra  in  terms  of  £'^Flux(E) 
in  order  to  highlight  the  energy  of  peak  power  output  in 
the  OSSE  band  and  emphasize  the  differences  between 
spectra  measured  at  different  flux  states.  The  spectrum 
with  open  squares  is  the  sum  of  all  days  with  a  flux  (in 


the  45-140  keV  band)  of  less  than  0.06  7cm“^  s“'.  The 
spectra  with  filled  squares,  open  circles  and  filled  circles 
are  for  a  flux  between  0.06  and  0.08,  between  0.08  and 
0.105  and  greater  than  0.105  7cm“^  s“',  respectively.  The 
low  amplitude  spectra  are  more  like  a  power  law  than  the 
other  data. 

5.  Search  for  line  features 

To  search  for  any  spectral  features  associated  with  specific 
flux  levels  that  would  have  been  averaged  out  in  the  to¬ 
tal  spectrum,  the  data  were  divided  into  individual  days 
and  each  day  was  analyzed  separately.  Following  the  anal¬ 
ysis  of  Ling  et  al.  (1987),  we  searched  for  the  presence  of 
MeV  bumps  by  fitting  each  day  of  data  with  a  thermal 
bremsstrahlung  spectrum  and  a  1.2  MeV  wide  gaussian 
line  feature  centered  at  1  MeV.  We  see  no  significant  am¬ 
plitude  for  a  broad  gaussian  feature  for  any  flux  level  of 
the  source.  When  all  the  data  are  added  together,  as  in 
Figure  1,  the  90%  confidence  level  (CL)  upper  limit  to  the 
1  MeV  line  is  4  x  10“^  7cm“^  s“',  compared  to  the  flux  of 
1.6  X  10“^  7cm“^  s“'  measured  by  Ling  et  al.  The  95%  CL 
upper  limit  for  a  1-1.5  MeV  wide  feature  near  4  MeV  is  2.0 
X  10“®  7cm“^  s“'  .  The  search  for  lines  can  be  repeated 


Fig.  3.  Photon  spectrum  multiplied  by  the  energy  squared. 
The  data  were  grouped  into  four  bins  according  to  their  45-140 


keV  flux. 


at  any  energy  and  width.  Given  the  previous  reports  of 
a  narrow  511  keV  line,  we  fit  the  data  for  a  narrow  (< 
instrument  resolution)  511  keV  line  for  each  day  of  data. 
We  do  not  find  any  days  with  a  significant  excess  511  keV 
emission  nor  a  correlation  between  the  511  keV  emission 
and  the  overall  source  strength.  The  90%  CL  upper  limits 
to  the  amplitude  are  between  0.5  and  1.0  x  10“^  7cm“^ 
s“'  for  most  days.  The  90%  CL  upper  limit  for  a  narrow 
511  keV  in  the  total  spectrum  is  7  x  10“®  7cm“^  s“'. 
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6.  Discussion 

The  variability  of  Cyg  X-1  at  gamma  ray  energies,  is  a  fac¬ 
tor  of  ~15,  much  greater  even  than  the  factor  of  ~3  mea¬ 
sured  during  the  HEAO-3  mission.  However,  the  variabil¬ 
ity  does  not  seem  to  be  between  discrete  states  but  rather 
among  a  continuous  range  of  possible  flux  values.  The  op¬ 
tically  thin  thermal  bremsstrahlung  fits  to  the  individual 
days  of  data  show  that  the  source  usually  fluctuates  in  a 
seemingly  random  fashion  in  a  narrow  range  of  intensity- 
temperature  phase-space.  This  is  unlike  the  behaviour  ob¬ 
served  in  GRO  J0422-1-32,  where  the  temperature  is  mono- 
tonically  increasing  in  time  as  the  transient  decreases  in 
brightness  (Kroeger  et  al.,  1995).  The  spectrum  is  best 
described  by  a  model  with  optically  thin  Comptonization. 
However,  a  preliminary  analysis  reveals  that  this  model 
cannot  be  applied  to  contemporaneous  GINGA  (Ebisawa, 
1994)  and  OSSE  data.  A  more  complex  model  using  an 
exponentially  truncated  power  law  with  an  angle  depen¬ 
dent  reflection  component  appears  to  fit  both  sets  of  data 
simultaneously  (Zdziarski,  1995).  Extrapolating  the  expo¬ 
nentially  truncated  power-law  model  to  lower  energies,  we 
find  that  the  calculated  flux  in  the  3-11  keV  range  for 
viewing  period  318  is  ~3  times  larger  than  for  a  typi¬ 
cal  observation.  This  anticorrelation  between  the  extrap¬ 
olated  low  energy  X-ray  flux  and  the  measured  gamma  ray 
flux  is  consistent  with  earlier  obervations  (Priedhorsky  et 
al.  1983)  and  is  not  very  model  dependent.  The  crossover 
energy  for  the  low  intensity  gamma  ray  spectra  to  become 
dominant  is  between  10-20  keV. 

We  have  searched  for  the  amplitude  of  a  broad  MeV 
bump  and  a  narrow  511  keV  line.  The  distribution  of  the 
daily  amplitudes  of  1  MeV  peaks  (figure  4)  is  gaussian, 
with  a  small  non-zero  mean,  indicating  the  underlying 
thermal  bremsstrahlung  distribution  underestimates  the 
continuum  flux  at  1  MeV  by  ~5  x  10“^  ph  cm“^  s“^. 
There  is  no  evidence  for  an  excess  of  points  at  high  flux 
values,  and  no  single  day  with  an  amplitude  near  that 
measured  by  Ling  et  al.  (1987)  over  a  12  day  period.  The 
90%  GL  upper  limit  to  the  1  MeV  bump  flux  is  in  the  1-4 
X  10“^  7cm“^  s“^  range  for  most  days.  If  MeV  bumps 
of  the  intensity  reported  by  Ling  do  exist,  they  must  be 
transient  within  the  7I  state.  Our  nondetection  of  MeV 
features  limits  the  duty  cycle  of  such  a  process.  For  1-day 
emission  episodes,  there  is  a  90%  probability  we  would 
have  detected  a  duty  cycle  greater  than  2.3%,  or  more 
than  25  outbursts  over  the  first  three  years  of  the  mis¬ 
sion.  For  4-day  episodes  the  duty  cycle  must  be  less  than 
4.5%,  and  for  14-day  episodes,  it  must  be  less  than  9%, 
which  corresponds  to  ~6.5  outbursts  over  the  mission  so 
far.  According  to  the  models  of  Dermer  et  al.  (1989)  and 
the  compilation  of  data  by  Bassani  et  al.  (1989),  there 
should  be  an  anticorrelation  between  the  hard  X-ray  flux 
and  the  gamma  ray  flux.  We  have  binned  our  data  in  two 
bands,  0.06-0.4  MeV  and  0.4-1. 5  MeV  for  a  direct  com¬ 
parison  with  Bassani  et  al.  From  the  total  data  set,  one 
would  conclude  that  there  is  a  correlation,  not  an  anticor¬ 
relation,  between  the  gamma  ray  flux  and  the  hard  X-ray 
flux,  consistent  with  higher  temperatures  at  higher  flux. 


-0.010  0.000  0.010  0.020 
Amplitude  (photons  cm'^  s'^) 

Fig.  4.  Distribution  of  the  daily  amplitudes  of  the  1  MeV  gaus- 
sians.  The  amplitude  reported  by  Ling  et  al.  (1987)  is  shown 
for  comparison. 

With  the  very  low  flux  observations  removed(45-140  keV 
flux  <  0.08  ph  cm“^s“^,  see  figure  2),  the  correlation  coef¬ 
ficient  of  this  ratio  is  0.03,  implying  no  correlation  between 
the  fluxes  during  most  observations. 
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